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We have studied the J/ip suppression in 158 GeV/c Pb+Pb collisions at CERN SPS. J/ip pro- 
duction is assumed to be a two step process, (i) formation of cc pair, which is accurately calculable 
in QCD and (ii) formation of J/tp meson from the cc pair, which can be conveniently parameter- 
ized. In a pA/AA collision, as the cc pair pass through the nuclear medium, it gains relative square 
momentum. As a result, some of the cc pairs can gain enough momentum to cross the threshold 
to become open charm meson, leading to suppression in pA/AA collisions. The model without any 
free parameter could describe the of NA50 data on centrality dependence of the ratio's; J/tp over 
I Drell-Yan, J ftp over minimum bias and also the Drell-Yan over minimum bias. The model was used 

to predict J/tp suppression at RHIC energy. At RHIC energy, hard processes may be important. 
With hard processes included, J/ip's are strongly suppressed, in agreement with other model calcu- 
lations. We also show that centrality dependence of J/tp over minimum bias ratio can be used to 
determine the fraction of hard processes in the collision. 
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I. INTRODUCTION 
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OO ■ One of the predictions of quantum chroniodynamics is the possible existence of a deconfined state of quarks and 
' gluons. In relativistic heavy ion collisions, under certain conditions (sufficiently high energy density and temperature) 
ordinary hadronic matter (where quarks and gluons are confined) can undergo a phase transition to deconfined matter, 
commonly known as quark gluon plasma. Over the years, nuclear physicists are trying to produce and detect this 
. new phase of matter, first at CERN SPS and now at Relativistic Heavy Ion Collider (RHIC). Indeed, it has even 
been claimed that evidence of new state of matter is already seen at Pb-|-Pb collisions at CERN SPS ||l|. The claim 
was based mainly on the experimental J/tp suppression and strangeness enhancement. Due to intrinsic ambiguities 
in disentangling similar effects arising from hadronic processes, the claim is considered to be rather controversial. 
O \ J/ip suppression is recognized as an important tool to identify the possible phase transition to quark-gluon plasma. 
5 ■ Because of the large mass of the charm quarks, cc pairs are produced on a short time scale. Their tight binding also 
makes them immune to final state interactions. Their evolution probes the state of matter in the early stage of the 
collisions. Matsui and Satz Q predicted that in presence of quark-gluon plasma (QGP), due to screening of color 
force, binding of cc pairs into J/ip meson will be hindered, leading to the so called J/ip suppression in heavy ion 
collisions . Over the years several groups have measured the J/ip yield in heavy ion collisions (for a review of the 



data and the interpretations see Refs. |Q,y). In brief, experimental data do show suppression. However this could be 
attributed to the conventional nuclear absorption, also present in pA collisions. 

The latest data obtained by the NA50 collaboration Q on J/ip production in Pb+Pb collisions at 158 A GeV is 
the first indication of the anomalous mechanism of charmonium suppression, which goes beyond the conventional 
suppression in a nuclear environment. The ratio of J/ip yield to that of Drell-Yan pairs decreases faster with Et 
in the most central collisions than in the less central ones. It has been suggested that the resulting pattern can be 
understood in a deconfinement scenario in terms of successive melting of charmonium bound states |5|] . Blaizot et al. 

have shown that the data can be understood as an effect of transverse energy fluctuations in central heavy ion 
collisions. Introducing a factor e — ET/Exib) , assuming that the suppression is 100% above a threshold density (a 
parameter in the model), and smearing the threshold density (at the expense of another parameter) the best fit to the 
data was obtained. Extending the Blaizot's model to include fluctuations in number of NN collisions at a fixed impact 
parameter, NA50 data could be fitted with a single parameter, the threshold density, above which all the J/tp mesons 
melt 0,11. Assumption that all the J/ip mesons melt above a threshold density, implicitly assume that a QGP like 
environment is produced in the collisions. NA50 data could also be explained in the conventional approach, without 
invoking QGP like scenario. Capella et al. Q analyzed the data in the comover approach. There also, the comover 
density has to be modified by the factor e. Introduction of this adhoc factor e can be justified in a model based on 
excited nucleons represented by strings |jlffl. Recently we have shown that the NA50 data is well explained in a QCD 
based model [|ll|, where formation of J ftp horn cc pair is suppressed, due to gain in the relative square momentum 
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of cc pair as it moves through a nuclear medium [|12| . Some of the cc pair can gain enough square momenta to cross 
the open charm threshold, thereby reducing the J /tp cross section. The model predicted large suppression of J at 
RHIC energy, in agreement with other model predictions. 

In the present paper we have extended our earlier analysis of NA50 data of transverse energy dependence of J/?/; 
to Drell-Yan ratio to include the transverse energy dependence of J over MB ratio and the DY over MB ratio. 
Recently Capella et al observed that on the average, a J/V' and also a DY event has Et ~ 3 GeV less than a 
minimum bias event. Et loss in a J/-!/)(DY) event, though small, can have larger effect beyond the knee of the Et 
distribution. We have also studied the effect of Et loss in J/i/; and in DY events. We confirm the analysis of Capella 
et al that in the ratio J/'0 over DY, Et loss effect gets cancelled. But the effect is manifested in the ratio of 
J /ip over MB as well as in the ratio of DY over MB. Next, we apply the model to predict Et dependence of J/V" 
suppression in Au+Au collisions at RHIC energy. At RHIC energy, hard scattering may be important. Suppression 
increases with hard scattering. It is also shown that J /if) over minimum bias ratio can be used to obtain fraction of 
hard scattering in the collision. 

The plan of the paper is as follows: in section 2, we briefly present the model of J /ip production and absorption in 
AA collisions. In section 3, we obtain the transverse energy distribution in Pb+Pb collisions. In section 4, relevant 
formula's for J /i), DY and MB cross sections will be noted. Analysis of NA50 data on the centrality dependence of 
the ratio's, J/i/'/MB, DY/MB and J/t/i/DY will be presented in section 5. In section 6 we will show predicted J/ip 
suppression pattern in Au+Au collisions at RHIC energy. Summary and conclusions will be given in section 7. 



II. MODEL OF J/i/) PRODUCTION AND ABSORPTION 

Qiu, Vary and Zhang |l^] proposed a model to describe the J/tp suppression in nucleon- nucleus/nucleus- nucleus 
collisions. For the sake of completeness, we will briefly describe the model. Qiu, Vary and Zhang assumed that the 
production of J /%p meson is a two step process, (i) production of cc pairs with relative momentum square g^, and (ii) 
formation of J mesons from the cc pairs. Step (i) can be accurately calculated in QCD. The second step, formation 
of J/V' mesons from initially compact cc pairs is non-perturbative. They used a parametric form for the step (ii), 
formation of J jip from cc pairs. The J/ip cross section in AB collisions, at center of mass energy -y/i was then written 
as, 

where runs over all parton flavors, and = + 4m^. The K factor takes into account the higher order 

corrections. The incoming parton momentum fractions are fixed by kinematics and are Xa — {\/ x\ + ^Q^ j s + xf)I'2. 
and Xb = {^^/xj^+4Q^/s — xf)/2. Quark annihilation and gluon fusion are the major sub processes for cc production. 
In the leading log, they are given by p^ ], 

a,,^c.(g^) = ^^(i + |)v/T^, (2) 

where as is the QCD running couphng constant and 7 = Am'l/Q'^. In Eq.^ i^ce— is the transition probability 
that a cc pair with relative momentum square evolve into a physical J/?/; meson. In a nucleon- nucleus/nucleus- 
nucleus collision, the produced cc pairs interact with nuclear medium before they exit. Observed anomalous nuclear 
enhancement of the momentum imbalance in dijet production led Qiu, Vary and Zhang [ pT| to argue that the inter- 
action of a cc pair with nuclear environment increases the square of the relative momentum between the cc pair. As 
a result, some of the cc pairs can gain enough relative square momentum to cross the threshold to become an open 
charm meson. Consequently, the cross sections for J /^p production are reduced in comparison with nucleon-nucleon 
cross section. If the J /ip meson travel a distance L, the transition probability Fcc->j/,p{q'^) in Eq.|l| will be changed 
to, 

Fcc->j/4,{q^) ^ F,s^j/i,iq^ + £^L), (4) 

with being the relative square momentum gain per unit length. Qiu, Vary and Zhang considered three different 
parametric forms (representing different physical processes) for the transition probability. All the three forms could 
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describe the experimental energy dependence of total J/tjj cross section in hadronic collisions ||Tl|. However, only the 
following form, 

could describe the experimental J lt\} data as a function of effective nuclear length jl^. In Fig.l, we have shown the fit 
obtained to NA50 experimental data on the J/^ total cross section as a function of effective nuclear length. We 
have used CTEQ5M parton distribution functions iQ. The parameter values, iirA^j/^=0.458, e^=0.225 GeV"^ / fm 
and aF=l are very close to the values obtained in Ref. 
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FIG. 1. Total J/ijj cross sections with the branching ratio to jJ.'^ fj,~ in proton-nucleus, proton-nucleus and nucleus-nucleus 
collisions, as a function of the effective nuclear length L(A,B). 



III. TRANSVERSE ENERGY DISTRIBUTIONS 



We have assumed that the probability to obtain Et at impact parameter b can be parameterized as, 

P(5, Et) cx exp{~{ET ~ qNj,{h)f /2q^aNp{h)) (6) 

The mean value of the transverse energy is < Et > (b) = qNp{h), where q is the average transverse energy per 
participant, and the dispersion is (t|;^ = aq^Np{b), where a is a dimensionless parameter. The number of participants 
Np{b) at impact parameter b is calculated as, 

Npib) = J dhnp{s,h) (7) 

where the transverse density np{h, s) is, 

n(b, s) = Ta(s)[1 - e-"""^«(*'-^)] + Teih - s)[l - e-"""^-^^")], (8) 

Nuclear thickness function, Ta,b{s) = / dzpA.si^, z) contain all the nuclear information. In the present calculation, 
we have used the following parametric form for 

PA{r) = , , ^"r-r, - (9) 

with a=0.53 fm, ro = I.IA^^^ . The central density is obtained from / pA{r)<Pr = A. 

The parameters q and a were obtained by fitting NA50 experimental transverse energy distribution. NA50 collab- 
oration B did not correct the Et spectra for the efficiency of target identification algorithm, which is lower than 



3 



unity for Et lower than 60 GeV. To obtain the parameters q and a we have fitted the purely inclusive part of the Et 
spectra {Et > 60 GeV). The fit is shown in Fig. 2. It was obtained with q=0.274 GeV and a=1.97. 
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FIG. 2. Experimental transverse energy distribution in Pb+Pb collisions along with the fit obtained assuming Gaussian 
Et — b correlation. 



IV. Et dependence OF J/t/;, DY AND MB EVENTS 



Here we present in brief, the relevant formula's for calculating minimum bias, Drell-Yan and J/ip cross sections. 
Minimum bias cross sections are easy to calculate. It is essentially the inelastic cross-section. In the Glauber model, 
at impact parameter b minimum bias cross section can be calculated as. 



d^'a"'" IdETd^h = (1 - exp{-aNNTAB{b))Pib, Et) 



(10) 



The Drell-Yan pairs do not suffer any final state interactions and the cross section at an impact parameter b as a 
function of Et can be written as, 



d^a'''' /dETd% = (7^]v j <fsTA{s)TB{s - b)P(6,ST), 



(11) 



where cr^^ is the Drell-Yan cross section in NN collisions. 

While Drell-Yan pairs do not suffer interactions with nuclear matter, the J/Tp mesons do. In the present model 
suppression factor depend on the length traversed by the cc mesons in nuclear medium. Consequently, we write the 
J l^p cross section at an impact parameter b as. 



d^a-^l'l' /dETd% = all^ j d'^sTA{s)TB{s - h)S{L{h,s))Pib, Et), 



(12) 



where ajj^ is the J/ip cross section in A'"A^ collisions and S{L(h, s)) is the suppression factor due to passage through 
a length L in nuclear environment. The length L{h, s) that the J/ip meson will traverse can be obtained as, 

L(b,s) = n(b,s)/2/9o (13) 

where n{b, s) is the transverse density [Eq.0|. Suppression factor ^(^(b, s) can be calculated using Eq.|l|, noting that a 
cc pair gains relative square momentum e , traversing a square of length L. Parametric value of e^, as shown before, 
was obtained by fitting nucleon-nucleus and nucleus-nucleus J/tp cross section data containing all Et- However, Eq.n2 
corresponds to a particular Et- Accordingly, square momentum gain factor needs to be modified. We modify the 
momentum gain factor to take into account the Et dependence as. 



e^ET)^e, 



,L{Et) 



(14) 
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where Eq is the momentum gain factor for all Et (which was obtained by fitting experimental data). L{Et) is the 
length through which a J/V' meson with transverse energy Et will travel. The length L{Et) can be calculated as [D, 

jd^hd\sTA{s)TB{h-s)[TA{s)+TB{h~s)]P{b,ET) 

^ 2poJdm^sTA{s)TB{h-s)P{b,ET) ■ ^ ' 

Fluctuations of transverse energy at a fixed impact parameter plays an important role in the explanation of the 
NA50 data. Above 100 GeV, i.e., approximately at the position of the knee, the 2nd drop in the data is due to the 
fluctuations in Et ■ In order to account for the fluctuations, following Capella et al we calculate, 

F{Et) = Et/E^^{Et). (16) 

where, 

jdHEr{b)P{b,ET) 

jd^hP{h^ET) ^'^^ 

The function F(Et) is unity up to the knee of the distribution, and increases thereafter, precisely where the 
fluctuations dominate. The replacement, 

L{h,s) ^ L{h,s)F{ET), (18) 
then properly accounts for the fluctuations in the Et distributions. 



V. COMPARISON WITH EXPERIMENTS 

Capella et al ||l^ observed that in a J/^p event sample and also in a DY event sample, Et ''^ 3 GeV is taken by 
the trigger and thus the transverse energy deposited in the calorimeter by the other hadrons is slightly less than the 
corresponding one in a minimum bias event sample. The Et loss in a J/^ event, though very small, can have larger 
effect in the ratio of J/ip over MB at the tail of the Et distribution. However, in the ratio of J/ip over DY the effect 
of Et loss cancels out. Exact magnitude of Et loss in a. J/ip event is difficult to estimate. Capella et al single out 
the NN collision in which J/ip (or DY) is produced and assumed that in that collision, no other hadrons are produced. 
We have implemented this in the present calculation by replacing the number of participants Np{b) in a J/^p or DY 
event by, 

Npib) Np{b) - 2, (19) 

minimum number of participants being two. 




Et (GeV) 

FIG. 3. Filled and square circles are the J/ip to Drell-Yan ratio in Pb+Pb collisions, obtained by NA50 collaboration in 1996 
and 1998 respectively. The solid line is the model prediction with or without Et loss effect. 
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In Fig. 3, we have compared the model prediction for the transverse energy dependence of J/i/" over Drell-Yan ratio, 
with the NA50 experimental data. We have assumed the normalization factor a'^/^ /a^^ =bi.h The solid line 
is the model prediction, obtained without Et loss consideration. Identical result is obtained when Et loss effect is 
taken into consideration. In the J /ip over DY ratio, the effect get cancelled to a high accuracy. The model gives 
good description to the NA50 data. The 2nd drop in the ratio at 100 GeV is reproduced. We note that without Et 
fluctuations, that drop is not reproduced. The result clearly shows that it is not essential to assume a deconfined 
scenario to explain the NA50 data on centrality dependence of J /ip to Drell-Yan ratio. Nuclear absorption alone is 
capable of explaining the data. 
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FIG. 4. NA50 data on the centrality dependence of the ratio DY/MB is compared with present model. The solid line is the 
model calculation without Et loss effect. The dashed line incorporate the Et loss effect. 



In Fig. 4, NA50 data on the centrality dependence of the Drell-Yan over minimum bias ratio is shown. The data 
were taken in 1996 using a thick target. NA50 collaboration does not consider the data reliable beyond the knee of 
the Et distribution, due to possible contamination of rescattering in the thick target. The solid line is a fit obtained 
to the data without Et loss effect. The ratio saturates beyond 100 GeV. When effect of Et loss in a DY event is 
taken into consideration (the dashed line) the ratio drops beyond 100 GeV. The effect is significant. At large Et, the 
ratio drops by more than 60% when Et loss effect is included. Indeed, experimental data do show such a drop, but 
as told earlier the data are not reliable beyond 100 GeV. The prediction can easily be verified in future experiments. 

NA50 collaboration reported the preliminary analysis of centrality dependence of J /^j) over MB ratio jl^. In Fig. 5, 
we have shown the preliminary data. The ratio increases up to the knee of the Et distribution, there after it decreases. 
There is an indication of change in slope around £'t=40 GeV. In Fig. 5, dashed and solid lines are model calculations, 
obtained with and without Et loss effect. The model predictions agree with experimental data. The change of slope 
around i?T=40 GeV is also reproduced. As with DY/MB ratio, the two predictions differ only beyond 100 GeV, where 
Et fluctuations dominate. With Et loss taken into account, the ratio drops faster, fitting the data better. Here also, 
the effect is significant, 60% or more at large Et- 

NA50 collaboration analyzed the J /ip data in two distinct manner. One is the standard analysis, in which the ratio 
of J I'll) over Drell-Yan (DY) was measured. In the other, which is called minimum bias (MB) analysis, the Drell-Yan 
events were replaced by the minimum bias events. The ratio of J /ip over minimum bias cross section was then divided 
by the ratio of theoretical DY over minimum bias cross section to obtain 
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The standard analysis data cover only up to the knee of the iJ^-distribution. The behavior of the J/ ip over DY 
ratio beyond the knee of the £'T-distribution is obtained entirely from the minimum bias analysis. Without Et loss 
effect, DY/MB remain practically constant beyond 100 GeV, and therefore the behavior of J/ip over DY and that 
of J/ip over MB is practically the same. However, present calculations clearly indicate that, due to Et loss effect, 
DY/MB does not remain constant but changes quite significantly in this Et range. Then behavior of J/ip over MB 
and that of J/ip over DY are not the same. True ratio will be flatter than obtained from the minimum bias analysis. 
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As suggested by Capella et al, it is better not to superimpose two analysis in a single figure. Models can be tested 
better by fitting them separately. 
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FIG. 5. NA50 (preliminary) data on tfie centrality dependence of tlie ratio J/?/)/MB is compared with present model. The 
solid line is the model calculation without Et loss effect. The dashed line incorporate the Et loss effect. 

VI. PREDICTION FOR RHIC ENERGY 

Present model can be used to predict centrality dependence of J/tp over DY in Au+Au collisions at RHIC energy. 
Experimental Et distribution at RHIC energy is not available, but we can make an educated guess for it. 
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FIG. 6. Simulated Et distribution for Au+Au collisions at RHIC energy is fitted with Gaussian Et — b correlation. 

Recent PHOBOS experiment showed that for central collisions, total multiplicity is larger by 70% at RHIC than 
at SPS. We assume that Et is correspondingly increased [Q. Accordingly, we rescale the Et distribution for Pb+Pb 
collisions and assume that it represent the experimental Et distribution for Au+Au collisions at RHIC (small mass 
difference between Au and Pb is neglected). At RHIC energy, the so-called hard component, which is proportional 
to the number of binary collision, appears. Model dependent calculations indicate that the hard component grows 
from 22% to 37% as the energy changes from y^=56 GeV to 130 GeV |2^. However we choose to ignore the hard 
component in i^T-distribution. The multiplicity distribution obtained by the PHOBOS collaboration, in the rapidity 
range 3 <| 77 |< 4.5 could be fitted well with or without this hard component. Indeed, it appears that the data are 
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fitted without the hard component |gO[]. Global distribution e.g. multiplicity or transverse energy distributions are 
not sensitive to the hard component. 

In fig. 6, filled circles represent the "experimental" Et distribution for Au+Au collisions at RHIC, obtained by 
scaling the Et distribution in Pb+Pb collisions at SPS. The solid line is a fit to the "experimental" i?T-distribution 
obtained with q—A6 GcV and a— 1.97. Nucleon-nucleon inelastic cross section (aNN) was assumed to be 41 mb at 
RHIC, instead of 32 mb at SPS 
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FIG. 7. Predicted centrality dependence of J/ip/UY ratio at RHIC energy for Au+Au collisions. The dashed and solid 
lines are present model predictions with and without hard scattering component in the transverse density. The dotted line is 
obtained in a model where all the J ftp's melts above a threshold density. 

While the Et distribution at RHIC energy may not depend on the so-called hard component which is proportional 
to number of binary collisions, J/ip suppression will strongly depend on the hard component, as it effectively increases 
the density of the nuclear medium. For / fraction of hard scattering, transverse density n(b, s) in Eq.^ is modified to 

@- 

n,w(b, s) ^ (1 - f)n{h, s) + fn'^'^^^h, s), (21) 

with ri'""''^(b, s) — (TNNTA{s)TB{h — s). With hard component, transverse density is increased, as a result, suppression 
will be increased at RHIC energy. In Fig. 7, solid line is the predicted centrality dependence of the J/tp/DY ratio 
without any hard scattering component in the transverse density. Suppression pattern is similar to that in Pb+Pb 
collisions at SPS energy. The 2nd drop now shifts at higher Et- With 37% hard scattering component in the 
transverse density (the dashed line) suppression increases and the effect of Et fiuctuations, i.e. the second drop is 
washed out. In Fig. 7, we have shown the prediction (dotted line) obtained in a model where all the J ftp's melts 
down above a threshold density The prediction also included 37% hard scattering. It agrees closely with the 

present model prediction (dashed line). Assumption of Jfip's melting above a threshold density implicitly assume 
QGP like environment is produced in the collision. Very similar prediction for Jftp suppression in a confined scenario 
and in a deconfined scenario, raises the possibility that it may not be possible to confirm the deconfinement phase 
transition, even at RHIC energy, from the Jfip data. Recently several authors have proposed that at RHIC energy, in 
a deconfined scenario, recombination of cc pairs will lead to enhancement of J/V''s, rather than its suppression [ p2[ . 
Inclusion of recombination effects may mask the large suppression obtained in a deconfined scenario. However, nuclear 
suppression as calculated presently will remain unaltered. It may then be possible to distinguish the deconfinement 
phase transition from the Jfip data. 

We have also calculated the centrality dependence of JftpfMB and of DYfMB ratio's at RHIC energy. Centrality 
dependence of DYfMB is very similar to that obtained in Pb+Pb collisions at SPS energy and is not shown. Centrality 
dependence of Jfip over MB ratio is shown in Fig. 8. The upper panel was obtained without any hard scattering 
component in the transverse density, while in the lower panel, hard scattering (37%) was included. As before, effect of 
Et loss is observed only beyond the knee of the Et distribution. It is interesting to note the difference in the shape of 
the ratio. Without hard scattering component, the ratio increases with Et up to the knee of the Et distribution, then 
decreases, a dependence similar to Pb+Pb data at SPS energy. With hard scattering included, the ratio increases at 
low Et^ then remains more or less same up to the knee of the Et distribution. Beyond the knee, it decreases again. 




8 



Very distinct shape of the two distributions can be used to determine experimentally the fraction of hard processes 
in the collision. 
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FIG. 8. Centrality dependence of J/?/>/MB at RHIC energy. Upper panel was obtained without hard scattering while 37% 
hard scattering is included in the lower panel. The solid and dashed lines are obtained without and with Et loss effect. 



VII. SUMMARY AND CONCLUSIONS 



We have studied the NA50 data on J ftp suppression in Pb+Pb collisions at SPS energy. The J/ip production is 
assumed to be a two step process, (i) production of cc pairs, which is accurately calculable in pQCD and (ii) formation 
of J/ip from the cc pair, essentially non-perturbativc, but can be conveniently parameterized. In the model, cc pairs 
gain relative square momentum as is travel through the nuclear environment, and some of the pairs can gain enough 
square momentum to cross the threshold to become an open charm meson. The parameters of the model were fixed 
by fitting experimental total J/iJj cross section in pp,pA and AA collisions. Without any free parameter, the model 
could very well describe the centrality dependence of the J/tp over DY ratio. The model also reproduces the centrality 
dependence of the J/ijj/MB and DY/MB ratio. We have also studied the effect of Et loss in a J/tp or DY event. 
Those events on the average produces ~ 3 GeV less Et than a minimum bias event. As shown by Capella et al [ p^ , 
the effect of Et loss get cancelled in the ratio J/tp/DY . But it affect the J/ip/MB ratio. Et loss also affect the 
DY/MB ratio. Ratio does not remain constant beyond the knee of the Et distribution, but drops. 

We have also used the model to predict the centrality dependence of the J /ip over DY ratio at RHIC energy. At 
RHIC energy, hard scattering component may be effective. J/-0 suppression depends strongly on the hard scattering 
component. Without any hard scattering, J /ip suppression closely follow the suppression pattern at SPS energy. The 
2nd drop occurring at higher Et- However, suppression increases when hard scattering is included and washes out 
the 2nd drop in the ratio. Suppression obtained in the nuclear environment closely agree with model prediction where 
QGP formation is implicitly assumed. Hard scattering also affect the J/i/j/MB ratio very strongly. While without 
any hard scattering, J /ip /MB increases up to knee of the Et distribution, with hard scattering included, the ratio 
remains more or less same over a large range of Et- The distinct change in shape raises the possibility of using 
J /ip/MB ratio to find the fraction of hard processes in the collision. 

To conclude, it is not essential to assume a deconfined scenario to explain the anomalous suppression of J /^ in 
Pb+Pb collisions at CERN SPS. The data are explained in terms normal nuclear suppression, also present in pa/AA 
collisions. Predicted suppression pattern at RHIC energy indicate, that even there, it may not be possible to identify 
deconfinement phase transition from the J /^ data. 

The author would like to thank Prof. A. Capella for his suggestions and comments. 
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